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Abstract 



I Multi-soliton form factors in sine-Gordon theory from the bootstrap are compared 

\^ ■ to finite volume matrix elements computed using the truncated conformal space ap- 

. proacli. We find convincing agreement, and resolve most of the issues raised in a 

previous work. 



1 Introduction 

The matrix elements of local operators (form factors) are central objects in quantum field 
theory. In two-dimensional integrable quantum field theory the S matrix can be obtained 
exactly in the framework of factorized scattering developed in [1] (for a later review see 
[2]). It was shown in [3j that in such theories using the scattering amplitudes as input it is 
possible to obtain a set of equations satisfied by the form factors. The complete system of 
form factor equations, which provides the basis for a programmatic approach (the so-called 
form factor bootstrap) was proposed in [4j . For a detailed and thorough exposition of the 
subject we refer to later this approach was also extended to form factors of boundary 
operators j6l [7] . 

Although the connection with the Lagrangian formulation of quantum field theory is 
rather indirect in the bootstrap approach, it is thought that the general solution of the 
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form factor axioms determines the complete local operator algebra of the theory. This 
expectation was confirmed in many cases by explicit comparison of the space of solutions 
to the spectrum of local operators as described by the ultraviolet limiting conformal field 
theory [8], [9l [TOl [HI [121 [El [HI [E]; the mathematical foundation is provided by the local 
commutativity theorem stating that operators specified by solutions of the form factor 
bootstrap are mutually local [5J . Another important piece of information comes from cor- 
relation functions. In the framework of quantum field theory, the operator matrix elements 
can be used to build a spectral representation for the correlation functions which provides 
a large distance expansion; this idea was implemented in integrable models using form fac- 
tors obtained from the bootstrap in [16]. On the other hand, the Lagrangian or perturbed 
conformal field theory formulation allows one to obtain a short-distance expansion, which 
can then be compared provided there is an overlap between their regimes of validity |17] . 
Other evidence for the correspondence between the field theory and the solutions of the 
form factor bootstrap results from evaluating sum rules like Zamolodchikov's c-theorem 
[l_8j or the A-theorem [19j, both of which can be used to express conformal data as spec- 
tral sums in terms of form factors. Direct comparisons with multi-particle matrix elements 
are not so readily available, except for perturbative or 1/A^ calculations in some simple 
cases [3j. 

In this paper we study form factors in finite volume, based on the approach developed in 
|20[[2T]. Finite volume form factors have also been studied in other approaches |22 | l23 | [M]: 
in addition, finite temperature form factors [25] are also related to this problem, as finite 
temperature is equivalent to compactified Euclidean time, and thus to a finite volume 
setting. 

One of the advantages of the framework developed in |20[ |2T] is that it allows for a 
direct comparison of solutions of the form factor axioms to field theory dynamics. This 
program has been successfully pursued in the case of diagonal scattering theories (those 
without particle mass degeneracies), both in the bulk and with boundary [2D1 [211 [2S1 127] . 
However, an extension to theories with non-diagonal scattering is still missing. The first 
steps were taken in [28] with a study of sine-Gordon breather and two-soliton form factors; 
later the effect of exponential corrections (more specifically so-called /i-terms) was also 
studied in detail l29j. Even earlier, finite volume breather form factors were already used 
in studying resonances [30] and form factor perturbation theory [31]. The present work 
is a natural continuation of this line of research, substantially extending and improving 
upon the previous results. Previously, there has been no way to study multi-soliton form 
factors because their integral representations could not be numerically evaluated. This was 
made possible by a quite involved and tedious numerical construction; the details of this 
technique are reported elsewhere |32) . 

It is important to realize that non-diagonal theories, whose spectra contain some non- 
trivial particle multiplets (typically organized into representations of some group sym- 
metry), such as sine-Gordon or the 0(3) nonlinear sigma model are very important for 
condensed matter applications (e.g. to spin chains; cf. [33]). The finite volume description 
of form factors can be used to develop a low-temperature and large- distance expansion for 
finite-temperature correlation functions [^[5^ 155]. which could in turn be used to explain 
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experimental data, e.g. from inelastic neutron scattering [5^ IM] . Another interesting 
application of finite volume form factors is the computation of one-point functions of bulk 
operators on a strip with integrable boundary conditions: the approach developed in |37] 
is in principle valid for general (i.e. non-diagonal) scattering. 

Therefore the extension to non-diagonal theories is an interesting direction. Sine- 
Gordon model can be considered as the prototype of a non-diagonal scattering theory, 
and it has the advantage that its finite volume spectra and form factors can be studied nu- 
merically using the truncated conformal space approach, originally developed by Yurov and 
Zamolodchikov for the scaling Lee- Yang model [38j, but later extended to the sine-Gordon 
theory [39]. Its exact form factors are also known in full generality |10t HH HJl HSl HI], and 
so it is a useful playground to test our theoretical ideas on finite volume form factors. 

To summarize, the motivations of the present work are: 

• To continue extending the description of finite volume form factors, initiated in [20| 
[21], to theories with non-diagonal scattering. 

• We perform the first detailed direct test of the sine-Gordon multi-soliton form fac- 
tors (conjectured from the bootstrap) along the lines of \'20\ m\ . i.e. by comparing 
them directly to numerically determined matrix elements computed from solving the 
explicit field theory dynamics in finite volume. 

• We also wish to make sure that the numerical representation developed in [32] are 
correct. These numerical results are intended to be used later for several independent 
lines of research, so testing and refining them is important. 

• In the previous work [28j, some issues were left unresolved. These were related to a 
sign observed in diagonal one-soliton matrix elements, and a numerical discrepancy 
in the comparison of diagonal matrix elements. Here we solve the first problem and 
present evidence that the second one is related to truncation errors inherent in the 
TCSA method. 

The paper is organized as follows. After a brief review of the necessary facts about sine- 
Gordon model in section |21 we recall the formalism for finite volume soliton form factors in 
section [3l Using the formalism developed in [28j, we give theoretical predictions for finite 
volume matrix elements between multi-soliton states in section [3l which are compared to 
numerical data from the truncated conformal space approach in section |H Section [5] is 
devoted to the conclusions and outlines remaining problems, to be investigated further. 

2 Brief review of sine-Gordon model 

2.1 Action and S matrix 

The classical action of the theory is 
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The fundamental excitations are a doublet of soliton/antisoliton of mass M. Their exact 
S matrix can be written as 111 



where 



s^io,0 = si^{e,OSo{e,0 (2.1) 



sxx{e,i) = s-_z{e,i) = i 
sXz{e,o = szt{0,O = ST{e,O 
s-1{9,o = s1-{9,o = Sr{9,o 



and 



sinh I 7 ) i sin ■ ^ 



sinh ( 1 sinh 



ITT- 



, . dt sinh 
J2sinh ^iizilie-"-?* + (e""-?* - l) (e-^-i)*/^ + e'^i^+Om) 



X 



2 sinh ^ cosh ^ 



(the latter representation is valid for any value of n G N and makes the integral repre- 
sentation converge faster and further away from the real 9 axis). Besides the solitons, the 
spectrum of theory contains also breathers; we omit details since these play no role in the 
sequel. We also introduced the parameter 



7r-/32 



Another representation of the theory is as a free massless boson conformal field theory 
(CFT) perturbed by a relevant operator. The Hamiltonian can be written as 

H = J dx^: (dt^)^ + {d^^f ■.+fxjdx: cos/3$ : (2.2) 

where the semicolon denotes normal ordering in terms of the modes of the fi = massless 
field. In this case, due to anomalous dimension of the normal ordered cosine operator, the 
coupling constant /i has dimension 

r 1 2-/32 Mtt 

jj, ~ [massj ' 

so it defines the mass scale of the model and the dimensionless coupling parameter is /3. 
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2.2 Soliton form factors 



We consider only exponentials of the bosonic field $. Their vacuum expectation value is 
known exactly |45] : 



Mv^r 



4lT 



2T 



/3V2 



"I 47r 


p dt 


exp < 


Jo ~ 



4n 



-2t 



sinh^ (i^t) 



2 sinh ( l^t ) cosh ( ( 1 - |^ ) t ) sinh t 



(2.3) 



with M denoting the soliton mass related to the coupling /i defined in via [IB] 

2-2A 

A 



2r(A) / v^r(^)M ' 

7rr(l-A)l 2r(^) 



Stt 



(2.4) 



As for multi-soliton form factors, at present there are three independent constructions: 
the earliest one by Smirnov (reviewed in [40j), the free field representation by Lukyanov 
|4H I4:2j and the work by Babujian et al. [ 43, .44j . Here we use formulae from Lukyanov's 
work |l2]; however, certain of his conventions are different and therefore we change the 
labeling of the form factors accordingly (see eqn. f l2.10p below). The reason is that the 
form factors we use satisfy form factor bootstrap relations which are slightly different from 
Lukyanov's conventions; in this we conform to the conventions of the papers [201 In 
our notations, the form factor equations are: 
I. Lorentz-invariance 



J^i + A, . . . , ^iv + A) = e^(^)^i^f....,(^i, . . . , ^iv) 



(2.5) 



where s{0) is the Lorentz spin of the operator O. 
II. Exchange: 



-^ii...ii,ife+i...ijv(^l' • • ■i(^k,Ok+l, . . . ,9n 



ll...«feife+l...«jV ^ 

^hh+r i^k - (^1, • • • , Ok,..., 9n) (2.6) 



III. Cyclic permutation: 

F?,.,...^^0, + 2z7r, ^2, . . . , e^) = e2-'-(°)i^J,.,^,^(02, • • • , ^i) (2.7) 

where oo{0) is the mutual locality index between the operator O and the asymptotic field 
that creates the solitons. 
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IV. Kinematical singularity 

Res + ^vr, e',e,,..., e^) = (2.8) 

e=e' 

where C is the charge conjugation matrix. 

V. Dynamical singularity 

- t Res + z«},/2 + e,e'- iu\J1, 6,,..., 6^) = rf,.Ff,^..,^ (9, 6,,..., 6^) (2.9) 

whenever k occurs as the bound state of the particles i and j, corresponding to a bound 
state pole of the S matrix, where T^j is the on-shell three-particle coupling and wf^ is the 
so-called fusion angle. The fusion angles satisfy 

ml = ml + + 2mimj cos u^j 
277 = u'lj + ui + u% 

and we also used the notation u^j = it — ufj. Equations I-V are supplemented by the 
assumption of maximum analyticity (i.e. that the form factors are meromorphic functions 
which only have the singularities prescribed by the axioms) and possible further conditions 
expressing properties of the particular operator whose form factors are sought. 
The form factors of the operator 

which satisfy equations (12. 5112. 8p with the locality index 

u}{Oa) = a mod 1 

can be obtained from 

-d...f72„(^l' • • • ?^2n) = (-l)''-^-a2iv -ai(^2n, • • • , 6*1) 

= (-ir47:U(^2.,...,^l) (2.10) 

where the functions T (derived by Lukyanov) are specified in appendix |X1 Equation (12. 9p 
for the dynamical singularities can then be used to construct form factors of breathers, 
which are bound states of a soliton with an antisoliton. 

Note that the sign factor (—1)" corresponds to a redefinition of the relative phase 
between a soliton and an antisoliton. In our previous work [28j, it was noticed that such 
a sign was necessary for a full agreement between the finite size data and the theoretical 
predictions. Since then we realized that this is explained by the difference between the 
conventions used for the form factor equations between Lukyanov's work jl2] and the 
finite volume form factor formalism developed in [201 1^ . Similarly, the other sign change 
is related to another difference in the conventions, namely the sign of the sine-Gordon field 
$, which can be compensated by either flipping the sign of a or exchanging the soliton 
with the antisoliton (charge conjugation). Finally, the rapidity ordering is again a matter 
of convention, this time that of fixing the basis for the asymptotic multi-particle states. 
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3 Soliton form factors in finite volume 



3.1 Finite volume form factors in non-diagonal theories 



The formulae for finite volume form factors, derived in pOl [2T] . were generalized for the 
case of non-diagonal theories in [28]. Here we only recall the necessary facts; for more 
details the reader is referred to the original papers. 

In finite volume L, the space of multi-soliton states can be labeled by momentum 
quantum numbers Ji, . . . , /jy. We introduce the following notation for them: 



IUl,/2 



(3.1) 



where the index r enumerates the eigenvectors of the n-particle transfer matrix, which can 
be written as 

r (A| {e„ e^})ii-x = '^Sf (a - ^i)5-^^^(a - e,) . ..sz^_^,^{x - ^^v) 

where 6i, . . . ,6n are particle rapidities. The transfer matrix can be diagonalized simulta- 
neously for all values of A: 



r (A| {6,, . . . , e^})^X ^n-,. = t^"-^ (A, {Ou}) ({^^}) 

We can assume that the wave function amplitudes \E''^^) are normalized and form a complete 
basis: 



«i...i]V 



these eigenfunctions describe the possible polarizations of the particle state with rapidi- 
ties $1, . . . ,6i\i inside the 2^ dimensional internal space indexed hj ii . . 

The rapidities of the particles in the state (13.11) can be determined by solving the 
quantization conditions 



Q 



Or, 



ML sinh 9j + 6^ 



27iL 



k 



(3.2) 



5^'\e^,...,9N) = -t\og&^ (%,R}) 



When considering rapidities which solve these equations with given quantum numbers 
/i, ... /at and a specific polarization state r, they will be written with a tilde as ^i, . . . , 6^. 
Using the above ingredients, the finite volume matrix elements can then be written as 



(^)({/(,...,/U|0(o,o)|{/i,...,M)f) 



p(r){9,,...,9^)p(^){9[, 



+ 0(e 



(3.3) 



where p^^^ and p^^^ denote the density of states of types r and s, 

ji---jM h---iN 

and the bar denotes the antiparticle. The absolute value in (13. 3p and in all similar formulae 
below is necessary to account for the different phase conventions of the multi-particle states 
used in the form factor bootstrap and in the finite volume calculations. 

Relation (13. 3 p is only valid for matrix elements with no disconnected pieces, i.e. when 
the rapidities in the two finite volume states are all different from each other. If there are 
particles with exactly coinciding rapidities in the two states, i.e. 9'j^ = 9i for some k and I, 
then there are further contributions. Note that equality of two quantum numbers such as 
I'f. = Ii is not sufficient for the presence a disconnected contribution, as the corresponding 
rapidities will in general be different due to the terms involving the phase shifts 6j^\ 
Therefore such terms are only present for the case when the two sets of quantum numbers 
are exactly identical, and also in the special case when the two states each contain a particle 
with exactly zero rapidity. At present, the disconnected terms are only known for states 
with diagonal scattering; the form of these contributions was obtained in |2T] . 

3.2 Soliton-antisoliton states 

This can be easily applied to soliton-antisoliton states. Two-soliton states form a four 
dimensional space corresponding to ss, ss, ss and ss. Due to the charge conjugation 
invariance of the S matrix, the eigenvectors of the transfer matrix have definite charge 
parity; together with charge conservation, this uniquely determines them. The eigenvectors 
of the two-soliton transfer matrix in the neutral subspace are [28] 

v^W = -^(0,+l,+l,0) 

vi/(-) = i=(0, +1,-1,0) (3.4) 

and are even/odd under charge conjugation, respectively. 

This results in the following quantization conditions for the soliton-antisoliton pair: 

Q\^\0i,92) = MLsinh^i + (5±(0i -^2) = 2nh 

Q\f\9i,92) = MLsinh^2 + ^±(^2 - ^1) = 27r/2 (3.5) 

where the phase-shifts 6± are defined from the eigenvalues of the two-particle ^'-matrix in 
the neutral subspace by 

S+{9) = SXZ{9) + S-l{9) = -e'^+(^) 
S_{9) = S+Z{9) - S-l{9) = e'^-^'^ 
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and the ± distinguishes the two states fl3.4p . Note the — sign introduced in the first hue 
which ensures that the phase-shifts are odd and continuous functions of the rapidity 9] as 
a consequence they vanish for 6* = 0. Due to this convention the + states are quantized 
with half-integer, while the — states are quantized with integer quantum numbers. 
The density of states can be written as the Jacobi determinant [28j 



From (13. 3p we obtain 



(0|O(0,0)|{Ji,/2}) 



(±) 

















802 




~9i - h 



(3.6) 



where 



in terms of (lA.ip and 61^2 are the solutions of (13. 5 p at the given volume L with quantum 
numbers Ii^2- This relation was already tested in [28j. 
Similarly one obtains 



W({/(,/^|O(0,0)|{/i,/2}) 



F^'\t7r + 9'2,tn + 9[,ei,e2] 



(r) 



p^^K0[,0',)p(r){9^,92) 



+ 0(e 



-fiL\ 



s, r = ±1 



(3.7) 



where 



F^'\9'„9[,9,,92y'^ 



9[,9i,92) 



'21 



F.++49',,9[,9,,92)+rF_+. 
sF+^+49'„ 9[,9i, 92) + rsF+__+(^^^, 9[,9r, 92 



provided the matrix element is non-diagonal, i.e. the two states differ either in their 
symmetry indices s,r, or in at least one of the momentum quantum numbers. 

States containing more than two solitons/antisolitons can be described using the al- 
gebraic Bethe Ansatz [28j; we do not enter into details as they are not needed in the 
sequel. 
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3.3 States containing only solitons 

Another way to test the multi-sohton form factors is to consider matrix elements where both 
states contain only solitons of like (say positive) topological charge. Since their scattering 
is diagonal, the formulae from |20[ |2T] are directly applicable. The quantization relations 
for these states are 



QuiOi, ...,9^) = MLsinhOk + - 9i) = 27ih 



k = l N 



(3.8) 



where the phase-shift S is defined by 

and the density of states is 
p(^i,...,ML = det J(^) 



, . . . , On] 



'kl 



09, 



kj 



N (3.9) 



Because of the sign in the definition of the phase-shift, the Ik are integer/half- integer for 
states containing on odd/even number of solitons, respectively. 

The finite volume matrix elements can be expressed as follows |20] : 

!({/(,..., /;,}|o(o,o)|{/i,...,j^})^l = 



fO(e-^^) (3.10) 

V(^^l,...,MLP(^i,---,^M)L 

Since the topological charge of the operator 

O = e^^* 

vanishes, the matrix elements are only nonzero when N = M. 

For the particular case of states containing only solitons we also know the form of 
disconnected contributions; since their scattering is diagonal, one can use the results from 
[21] . For diagonal matrix elements 

{{h,...,lN}\O{0,0)\{h,...,lN})L ^ 



P({i,...,iv})I^ ^'-''^ 

J2 -F(A)iP({l,...,iV}\A)i + 0(e-^^) 

Ac{l,2,...N} 

where 1^41 denotes the cardinal number (number of elements) of the set A 

p{{ki, kr})L = p{Ok,, • • • , Okr)L 
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is the r-particle Bethe-Yang Jacobi determinant ( 13. 9 p involving only the r-element subset 
1 < ki <■■■< kr < N of the N particles, and 



J'{{ki,...,kr})L 

F^{9i, . . . , 9i)i 



Fri^ki, ■ ■ ■ , dkr) 

lim + ivr + e, . . . , 01 + ZTT + e, 01, . . . , Oi) 

e-*-0 



Besides diagonal matrix elements, the only other possibility for disconnected terms to occur 
is when both states contain a stationary particle; in our case it can only happen in matrix 
elements with the same (odd) number of solitons on both sides. The general formula can 
be found in [21]; here we only quote the case needed in the sequel: 



l({r,o,-r}|o|{/,o, -/}),! = 

J-i,i(0'|0) + ML F^_^^{in + e\ ITT - 6', -e, 6) 



(3.12) 



+ 0(e 



-fiL\ 



where 



J-i,i(0'|0) = liniFf_ 



p(0',O,-0Olp(0,O,-0)l 

' + e,iiT -6' + e,iTT + e,0, -6, 



[ITT 



4 Numerical results 
4.1 Numerical methods 

To evaluate the form factors numerically, we use the truncated conformal space approach 
(TCSA) pioneered by Yurov and Zamolodchikov [38j. The extension to the sine-Gordon 
model was developed in [39j and has found numerous applications since then. The Hilbert 
space can be split by the eigenvalues of the topological charge Q (or winding number) and 
the spatial momentum P, where the eigenvalues of the latter are of the form 

27rs 

s is called the 'conformal spin'. The basis of the Hilbert space is constructed in the ultra- 
violet limiting massless free boson CFT with central charge c = 1, and a (dimensionless) 
upper cutoff E^ut is imposed on the scaling dimension (which is the sum of the left and 
right conformal dimensions) of the states kept under the truncation. 

In sectors with vanishing topological charge, we can make use of the symmetry of the 
Hamiltonian under 

C : <l>(x,t) ^ -$(x,t) 

which is equivalent to conjugation of the solitonic charge. The truncated space can be split 
into C-even and C-odd subspaces that have roughly equal dimensions [28j , which speeds up 



11 



the diagonalization of the Hamiltonian by roughly a factor of eight (the required machine 
time scales approximately with the third power of matrix size). We used the program 
developed for the work [28], with cutoff values Ecut ranging between 16 to 26; the highest 
cutoff was chosen such that the dimension does not exceed 11000 states (in order for the 
program to fit into available computer memory and also finish in a reasonable amount of 
time); the maximum E^ut permitted by this criterion depends on the value of the sine- 
Gordon coupling P, and the topological charge and spin of the sector under consideration. 
For the matrix element calculations, we chose the operator 

O =: e^'^* : 

which is essentially one half of the interaction term in the Hamiltonian in f l2.2p . The 
semicolons denote normal ordering with respect to the A = free massless boson modes. 
This operator has conformal dimension 

Ao = Ao = 

8n 

Using relation (12. 4p we can express all energy levels and matrix elements in units of (ap- 
propriate powers of) the soliton mass M, and we also introduce the dimensionless volume 
variable / = ML. The general procedure is the same as in [201 121] : the particle content 
of energy levels can be identified by matching the numerical TCSA spectrum against the 
predictions of the Bethe-Yang equations. 

To generate the data used for comparison, altogether 342 TCSA Hamiltonians were 
diagonalized, and from them 1350 operator matrices were computed. Three values of 
couplings were used: ^ = 2/7, 50/239 and 50/311; for each of them we evaluated the 
sectors Q = 0,1,2,3 with spins s = 0,1,2 (for the sectors with Q = 0, this was done 
separately for the C-even/odd projections) and with as many values for Ecut > 16 as the 
dimensionality constraint admitted. This left us with a vast amount of useful data of which 
we only include an illustrative sample; we performed the comparison for a much larger set, 
with identical results to the ones presented below. 

As in all our previous works on finite volume form factors (see e.g. [20| |2T[). energy 
levels predicted by the Bethe-Yang equation were used to identify the particle contents 
of the finite volume energy levels computed numerically from the TCSA method. Due 
to level crossings, at certain values of the volume L there can be more than one TCSA 
candidate levels for a given Bethe-Yang solution; for the data presented here we kept only 
unambiguously identified levels. 

In all of the figures presented below, we denote the operator matrix element by /: this 
means taking the absolute value of the matrix elements which is normalized by choos- 
ing the TCSA vectors orthonormal. This conforms to the conventions used in eqns. 
(I3.3|3.10|3.11|3.12p . In all cases, the discrete points are the numerical TCSA data, while 
the continuous lines are the corresponding theoretical expectations. 
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4.2 Sources of deviations 



There are two sources of deviations between the theoretical predictions and numerical 
results: 

1. Exponential finite size effects are neglected in the theoretical description for the vol- 
ume dependence of matrix elements, outlined in section [3l While they are partially 
understood (especially the so-called /i-terms [17 t 1 ^ ). there is no systematic descrip- 
tion for them yet, so we do not consider them here. Generally, they are expected to 
be larger for smaller ^: both because the breathers become lighter in terms of the 
mass scale M provided by the soliton mass (which affects so-called F-terms arising 
from breather loops non-trivially wound around the finite volume), and also because 
they become less tightly bound when considered as bound states of other breathers 
(which enhances the /i-terms related to compositeness) . Some of the /i-terms can 
also be dangerously enhanced by the analytic behaviour of form factors [29], but no 
sign of such behaviour was seen for the matrix element considered in this work. 

2. Truncation errors introduced by TCSA, on the other hand, generally increase with 
the volume and are also larger for higher excited states. In sine-Gordon theory, they 
have been observed to become smaller when decreasing ^, so the two sources of devi- 
ations behave the opposite way when the sine-Gordon coupling is varied. Behaviour 
of truncation errors in the asymptotic regime of large values of the cutoff can be 
theoretically described by a Wilsonian renormalization group |48[ HHl I5U] . 

Level crossings also present a problem in numerical stability, since in their vicinity 
the state of interest is nearly degenerate to another one. Since the truncation effect 
can be considered as an additional perturbing operator, the level crossings are even- 
tually lifted. However, such a near-degeneracy greatly magnifies truncation effects 
on the eigenvectors and therefore the matrix elements [27j. This is the reason behind 
the fact that there are some individual numerical points that are clearly scattered 
away from their expected place (cf. fig |4.ip . 

For any quantity (energy levels and matrix elements) which is compared between the 
theoretical predictions and the numerics, one can define the "scaling regime", which is the 
volume range in which the two sources of deviations are the smallest, i.e. the range in which 
truncation errors and exponential finite size effects are approximately the same magnitude. 
This range depends on the following factors: 

1. The value of the sine-Gordon coupling /3: when /3 (or equivalently) ^ decreases, it 
shifts to larger values of the volumes, and also becomes longer. 

2. The TCSA truncation: it becomes longer when increasing the value and also shifts 
to slightly higher values of the volume. 

3. The quantity under considerations: as shown below, diagonal matrix elements are the 
ones most affected by truncation errors, for which we have no theoretical explanation 
at present. 
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Figure 4.1: Non-diagonal form factors in the Q = sector for ^ = 50/239. A few examples 
of individual data points affected by truncation effects magnified by the vicinity of a level 
crossing can be seen at ML = 23, cf. the second and fourth lines from below (plotted with 
diamonds and stars, respectively). 

4.3 Results for four-particle form factors 

We can consider off-diagonal matrix elements between soliton-antisoliton two-particle states 
(for the diagonal ones we do not have the theoretical description yet, cf. the discussion in 
the conclusions). The theoretical prediction is given by eqn. (13. 7p and the comparison is 
shown in figure 14.11 

In addition we can use the soliton-soliton two-particle states. For off-diagonal matrix 
elements it is straightforward to use (I3.10p and we obtained a good agreement as demon- 
strated in fig. 14.21 For the diagonal case however, one observes a discrepancy between the 
predictions from eqn. (I3.1ip and the numerical results in fig. 14.31 which becomes smaller 
for smaller values of ^. As illustrated in fig. 14. 4^ this can be explained by truncation errors, 
which are indeed improved by decreasing ^. One can try to extrapolate the truncation de- 
pendence; however, it turns out that it does not fit the theoretically expected asymptotics 
derived in [19], which means that the leading order renormalization group behaviour is not 
yet valid at the cutoffs considered. Extrapolations reproducing the theoretical predictions 
can be found, but for a cut-off dependence which has an exponent that differs from the 
predictions of the renormalization group; in addition, the available range of cutoff values 
is not sufficient for a reliable determination of the exponent from the numerical data. In 
the conclusions we discuss how the situation can be improved, but this is out of the scope 
of the present work. 
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Figure 4.2: Non-diagonal form factors in the Q = 2 sector for = 50/239 



4.4 Six-soliton form factors 



We tested six-soliton form factors by comparing the predictions from eqn. (I3.10p to off- 
diagonal matrix elements between states composed of three solitons (and no antisolitons). 
The agreement is again very convincing, as demonstrated in fig. 14.51 For diagonal matrix 
elements, we noticed similar discrepancies as in the case of diagonal four-soliton form 
factors; however, the truncation dependence proved to be much worse in this case, so while 
the results were qualitatively consistent, they were not as good as for the four-soliton case. 

In addition, for this case there is an interesting new possibility of having disconnected 
parts originating from particles with exactly zero rapidity, described by eqn. (I3.12p . For 
these matrix elements we get a very convincing agreement again, as demonstrated in fig. 

mi 



5 Conclusions and outlook 

In this work we compared the conjectured exact soliton form factors of sine-Gordon theory, 
obtained from the bootstrap, to finite volume matrix elements given by the truncated 
conformal space approach. 

For non-diagonal matrix elements we find excellent agreement between the numerical 
results and theoretical expectations, both for four-soliton and six-soliton form factors. For 
diagonal matrix elements we found some discrepancy similar to the one noticed in [28] 
for four-breather form factors. This discrepancy tends to decrease for smaller value of ^ 
(or equivalently /3) and we argued that it can be attributed to truncation effects. Similar 
effects were observed for boundary form factors in pT] and based on the accumulated data 
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Figure 4.4: Truncation dependence of the diagonal matrix element 
({1/2, -l/2}|C|{l/2, -1/2}) in the Q = 2 sector for ^ = 2/7 
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Figure 4.6: Form factors with zero-momentum disconnected pieces in the Q = 3 sector 



we are inchned to think that there is something special about cutoff dependence of diagonal 
matrix elements. 

Convergence of the TCSA can be improved by renormalization group methods |18| 
Sni EDJ; as we discussed in subsection 14.31 it turns out that the TCSA data are not yet 
in the regime where the leading RG behaviour is applicable, and the extrapolation fits 
are not reliable enough to determine the exponent of the cutoff dependence. This can be 
helped by applying the numerical RG technique proposed in [19]; however, developing a 
systematic program for that takes a substantial amount of effort and time, and work in this 
direction has just started. One can also extend the domain of comparison by improving 
the theoretical description for smaller volumes (where truncation errors are negligible) by 
describing exponential finite size effects. 

Aside from the above-mentioned technical issues, there is a crucial missing piece, namely 
the description of disconnected pieces for states in which the scattering is non-diagonal, 
i.e. an extension of formulae ( 13.11^ I3.12p to the general case. The work aimed at resolving 
this issue is in progress, and the developments in this paper are also useful in preparing a 
testing ground for future theoretical conjectures. Once this final piece is in place, it will 
be possible to use the systematic formalism developed in [35] for the form factor expansion 
of finite temperature correlators to evaluate correlators in field theories with non-diagonal 
scattering, such as sine-Gordon theory or the 0(3) nonlinear cr-model. 
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A Explicit formulae for the soliton form factors 

Let us denote the form factor functions defined by Lukyanov by 
where ai = ± and 

2n 

fc=i 

is necessary for the matrix element to be different from zero. The operators Z are given 
by 



Z_{e) = [ ^^{i-2a-s./p-)[^-e)^-i-m^im 

V 4Ci 1^ Jc^ 27r 

-e ^ -^i 
Jc. 27r 



_g(l-2a-87r//32)(7-e)gi0(e)g-i<^(7) 



The averages are computed by the multiphcative Wick theorem (vahd for exponential 
operators) using 

The function G is given by the integral representation 

V2y 170 ^ VV y sinh2tcoshtsinht{J 

= iC, sinh (^0 f[g{9, C, kf exp | ^e"^^* (l + TV - TVe-^*) 



xsinhMfl-^^t^ -ht(e-l) 



Try y sinh 2i cosh i sinh 



19 



where the second formula provides an extension for the domain of convergence of the 
integral by factorizing out the following pole factors 



2 / , ,^ .„\ , (2fc~2+g)7r+i6) 



2 



p ^2fe-K j p ^ (2fc+g)7r-ie j p ^( 

p / (2fc-l)7r+t6> \ p ( 2fc-l+g \^ p (2fc-l+g)7r+ie \ 



r (f ) r (Mi^) r 



and is independent of A^, 



r_ prft sinh^(t/2)sinh(t(e-l)) l 
Ci _ G( Z7rj-exp| t sinh(2t) cosh(t) sinh(tO J 



C2 = exp<j4 

and 



' rftsinh^(t/2) sinh(t(^- 1)) 
~t sinh(2t) sinh(i^) 

1 



jY Y ( (2fc-5/2+7r)g+ie \ p f {2k~l/2)-n-ie \ p f 2fc-l/2 

2 TT V ^€ / V ^€ / V € 

cosh p ^-j^ _|_ 2fc-3/2 j ^ p ^ (2fc+V2)7r-ie j p ^ (2A;-3/2)7r+t6) ^ 



00 



X exp ^ -2 / ^e-^* sinh^ f f 1 - - '"^ ' " 



TT / / sinh 2t sinh 



sinh — sinh ^ 

4^ e 



where, again, the integral formula is eventually independent of the natural number A^; 
it provides a representation which converges faster numerically and is valid further away 
from the real Q axis with increasing N . The contours in the integrals are such that the 
"principal poles" of the ly-functions are always between the contour and the real line, 
where the "principal pole" of Wix) is the one located at s = —in/ 2. 

The integral representation can be evaluated in a closed form at the free fermion point 
,^ = 1 and also for the two-particle case when a is either integer or half-integer ^42j. Here 
we only quote the case needed in the text: 

G{e) /ttA 42cosh(f)e^^ 



:FU{e) = cot ^ (a.i) 

The numerical evaluation of the integral representation is rather involved; the details are 
given in [32]. 
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